Hidden vortices in a Bose-Einstein condensate in a rotating double-well potential 
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We study vortex formation in a Bose-Einstein condensate in a rotating double-well potential. 
Besides the ordinary quantized vortices and elusive ghost vortices, "hidden" vortices are found 
distributing along the central barrier. These hidden vortices are invisible like ghost vortex but carry 
angular momentum. Moreover, their core size is not given by the healing length, but is strongly 
influenced by the external potential. We find that the Feynman's rule can be well satisfied only 
after including the hidden vortices. There is no critical rotating frequency for the formation of 
hidden vortex while there is one for the formation of ordinary visible vortices. Hidden vortices can 
be revealed in the free expansion of the Bose-Einstein condensates. In addition, the hidden vortices 
in a Bose-Einstein condensate can appear in other external potentials, such as a rotating anisotropic 
toroidal trap. 

PACS numbers: 03.75.Lm, 03.75.Kk, 67.85. De 
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I. INTRODUCTION 

Double- well (DW) potential is an important model po- 
tential for its simplicity and yet richness in physics. The 
properties of Bose-Einstein condensates (BECs) in a DW 
potential have been studied extensively since the realiza- 
tion of BECs [1] . With the advances of technology, DW 
potentials for ultracold atoms can now be realized in ex- 
periments with great controllability and precision 043- 
One particular interesting development is the possibility 
to rotate a DW potential via radiofrequency dressing [5i] . 
This opens the door to study the behavior of vortices of 
degenerate quantum gas in a rotating DW potential. As 
topological defect, quantized vortices have contributed 
greatly to reveal the phase coherence, superfluidity and 
nonlinear phenomena for degenerate quantum gases, and 
have been the subject of extensive experimental and the- 
oretical studies [JH^. 

The DW potential also offers a unique testing ground 
for the well-known Feynman's rule of vortex [23]. The 
Feynman's rule is very powerful relation that links the 
total number of vortices with rotation angular frequency. 
Feynman argued that for a rotating superfluid with an- 
gular frequency fl, the superfluid should be regarded as 
a classical fluid when it reaches the steady state. This 
leads to an important mathematical relation that the to- 
tal number of vortices Ny in an area A is linearly propor- 
tional to fl, 27rhNy/m = 2^A. Alternatively, the Feyn- 
man's rule can expressed as l^/h — Ny/2 with being 
the mean angular momentum per atom at equilibrium 
[3, m, [2^, [2J] . This rule was originally for a uniform 
superfluid heliurn, and has been intensively studied both 
theoretically [Mill and experimentally %\X^ for a BEC 
trapped in single harmonic potential. It is interesting to 
know how the rule fares in a more complicated geometric 
confinement. The DW potential provides a clear oppor- 



tunity to answer this question. 

In this paper we have conducted a comprehensive two- 
dimensional (2D) numerical study of vortex formation 
in a BEC in a rotating DW potential. We find sur- 
prisingly that the in situ density distribution seems to 
violate significantly the Feynman's rule in that the to- 
tal number of vortices with visible core is significantly 
smaller than 2lz/h. With the belief that the Feynman's 
rule should hold in some form, we carefully analyze the 
results and find these "missing" vortices are distributed 
along the central barrier of the DW potential. Unlike the 
usual vortex, these vortices have no visible cores but have 
phase singularities, and their core size is not given by the 
healing length but is strongly influenced by the external 
potential. For this reason, they may be called "hidden 
vortices" . With the inclusion of the hidden vortices, one 
recovers the Feynman's rule. 

These hidden vortices remind us of ghost vortices found 
numerically in Ref . f23| . The ghost vortices always lie at 
the outskirts of the condensate where the particle density 
■01^ is very small. As a result, similar to hidden vortex, 
they show up in numerical results as phase singularities 
but have no visible cores. However, there are key dif- 
ferences: the ghost vortex carries no angular momentum 
while hidden vortex does; the core size of ghost vortex is 
determined by healing length like the usual visible vor- 
tex while that of hidden vortex by the shape of external 
potential. In addition, we find that with the increasing 
of rotation frequency f2 the hidden vortices appear first 
in the DW system, followed by ghost vortices and usual 
visible vortices. Furthermore, the angular momentum 
can be put gradually into the BEC via the generation of 
the hidden vortices while the emergence of a visible vor- 
tex is still accompanied by a jump in the system angular 
momentum. Although the hidden vortices are invisible 
in the in situ density distribution, after free expansion 



2 



of the BEC, they can appear in the density distribution 
because of their stable topological structure. 

It is well known that there exists a special type of 
vortex called Josephson vortex (or fluxon) in a long su- 
perconducting Josephson junction (26r] or between two 
weakly-coupled BECs 13, 2^ ■ These Josephson vortices 
can be regarded as hidden vortices. However, hidden vor- 
tex is a more general notion than the Josephson vortex 
(fluxon) as hidden vortex can exist in non-DW poten- 
tial (or non- Josephson structure). We have used rotating 
anisotropic toroidal trap to illustrate this point. 

This paper is organized as follows. In Sec. II, we 
present a phenomenological model to describe the dy- 
namics of a BEC confined in a rotating DW potential in 
the presence of dissipation. Hidden vortices are found in 
the rotating DW BEC, where the Feynman's rule is satis- 
fied only after including these hidden vortices. A simple 
and feasible scheme is proposed to observe the hidden 
vortices. In Sec. HI, we study the vortex formation pro- 
cess and the critical rotating frequency in the rotating 
DW BEC. In Sec. IV, we discuss the hidden vortices in 
a BEC confined in a rotating anisotropic toroidal trap. 
Sec. V provides a summary and discussion. 



II. HIDDEN VORTICES IN A ROTATING DW 
BEC AND FEYNMAN'S RULE 

We consider the situation where the condensate is 
tightly confined in the axial direction {uix,i^y <^ uJz) so 
that the system is effectively two dimensional. The DW 
potential is described by 



Vdw{x,v) 
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where Vq and a are the height and width of the poten- 
tial barrier, respectively, and A = Ljy/u)x denotes the 
anisotropy parameter of the harmonic trap. In the pres- 
ence of dissipation, the order parameter in the frame 
rotating with the angular velocity O around the z axis 
obeys the time-dependent Gross-Pitaevskii (CP) equa- 
tion 



-{'^l + ^l) + Vd^ + c\i;\' -nL, V- (2) 



Here — i{ydx — ^dy) is the z component of the 
angular- momentum operator, 7 characterizes the degree 
of dissipation, and c is the 2D interatomic interaction 
strength. In this work, length, time, energy, angular mo- 
mentum, and rotation angular frequency are in units of 
do = y^h/2mLUx, i/^x, fi^x, ft, and ujx, respectively. The 
phenomenological dissipation model ([2]) is a variation of 
that in Ref. [24] . For the case of a BEC in a rotating har- 
monic trap, our computation results agree well with the 
experimental observations in Ref. Q and the simulation 
results in Ref. 24]. 

In our calculations, we first obtain the initial ground- 
state order parameter in the DW potential by imaginary 



time propagation method [29''-'3l|] for J7 = 0. The vortex 
formation process is then studied by solving numerically 
Eq. ^ with different ft. Here we consider the BECs 
of *^Rb atoms with repulsive interaction. The system 
parameters are chosen to be = Wj, = 27r x 40 Hz, 
^ 2tt X 800 Hz, Vo = 40, ct = 0.5, c = 600. In 
Eq. ([2]), the variation of 7 only infiuences the relaxation 
time scale but does not change the dynamics of vortex 
formation and the ultimate steady vortex structure. In 
our computation, we choose 7 = 0.03, which corresponds 
to a temperature of about O.lTc (3^ . 

Fig. 1(a) shows the steady density distribution jV'l^ 
a.t t = 250 for the DW potential rotating with = 0.9. 
From this in situ density distribution, we see a pair of 
ordinary vortex lattices with triangular structure as ex- 
pected from the rotating DW configuration. However, 
by simply counting, we find something surprising. In 
Fig. 1(a), the total number of vortices is iV„ — 18. Nu- 
merical results show that the mean angular momentum 
per atom 1^ — J J ip* Lz4'dxdy / J J dxdy is about 
Ij, « 16 3> Ny/2. It seems that the Feynman's rule 
(3 [isl [23I , [2^ is no longer satisfied. For other rota- 
tion frequencies, this seemingly significant violation of 
Feynman's rule is also found. With the belief that the 
Feynman's rule should always hold, this violation indi- 
cates that some angular momentum is missing and is not 
manifested in the form of ordinary vortices. 

To find the missing angular momentum, we look into 
the phase distribution of ijj^x.y.t ~ 250), which is plot- 
ted in Fig. 1(b). We find that besides the phase singu- 
larities corresponding to the above mentioned vortices, 
there are other phase defects, distributed along the cen- 
tral barrier and the outskirts of the cloud. The initial 
reaction is that these phase defects, which are invisible 
in the in situ density distribution, are ghost vortices dis- 
cussed in Ref. ^24] . Since ghost vortices are known not 
carrying angular momentum, it seems that these invisi- 
ble phase defects can not account for the missing angular 
momentum. However, a more careful examination shows 
otherwise. 

The phase singularities along the central barrier are 
not ghost vortices and they carry angular momentum. To 
see this, we assume Nh is the total number of phase sin- 
gularities along the central barrier, and Nt = {Ny + Nh) 
is the sum of iV„ and Nh- From Fig. 1(b), we have 
Nh — 16. If we include them, the Feynman's rule 
Iz — Nt/2 — {Ny + Nh)/2 is well satisfied, indicating that 
these phase singularities do carry angular momentum like 
the usual visible vortices. We have checked other rotat- 
ing frequencies and that the Feynman's rule can always 
be well satisfied by including the phase defects along the 
central barrier. In Fig. 1(c), we have plotted the depen- 
dence of Iz on Nt/2 for different rotation angular frequen- 
cies. The solid line is the Feynman's rule Iz — Nt/2 while 
the solid circles are the numerical results of Iz at equi- 
librium. The excellent agreement between them strongly 
support that phase singularities along the central barrier 
carry angular momentum, and thus are not ghost vor- 
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FIG. 1; (color online) (a) Density distribution \Tp\'^ and (b) 
phase distribution of ip a.t t = 250 after rotating the system 
with Q — 0.9. (c) Iz versus Nt/2 for different Jl, where the line 
is Iz — Nt/2 = {N-u + Nh) /2 and the solid circles denote the 
numerical results. The value of the phase varies continuously 
from to 27r. The darker color area indicates the lower density 
or phase. The small derivation between the line and the circles 
in (c) may be due to the inhomogeneous condensate density 




tices. We call these topological defects along the central 
barrier "hidden vortices" for their difference from the or- 
dinary vortex and the ghost vortex. 

Two factors are involved in why a hidden vortex car- 
ries angular momentum while a ghost vortex does not: 
location and core size. To see this, we consider a phase 
defect with a singly quantized circulation in the conden- 
sate. The angular momentum carried by this phase defect 
varies with the location as Iz ~ [1 — / B?) where 
R is the size of the condensate and r is the distance from 
the center. Since a ghost vortex always lies on the out- 
skirts of the condensate, meaning r ~ i?, its contribution 
to angular momentum is negligible. For a hidden vor- 
tex, which locates near the center of the condensate, we 
have r < R. Therefore, its contribution to the angular 
momentum is significant and needs to be counted for the 
Feynman's rule. Furthermore, the core size of a non- 
hidden vortex is about the healing length ^ = -y/l/nc (n 
is the local density of the condensate without vortices). 
Since n is very small for a ghost vortex, the core size of 
a ghost vortex approaches an infinite value. As a result, 
ghost vortices neither contribute to the angular momen- 
tum nor the energy of the system. For a hidden vortex, n 
is also very small at its location. However, its core size is 
determined by the barrier width, not the healing length. 
Therefore, hidden vortex can contribute to the angular 



momentum. With the local density approximation, our 
numerical calculations do show that the hidden vortices 
carry significant angular momentum while the angular 
momentum due to a ghost vortex can be negligible. 

Even though the hidden vortices are invisible in the in 
situ density distribution as shown in Fig. 1(a), we find 
numerically that they show up in the cloud after free 
expansion (see the following discussion). This makes it 
possible to observe and test the existence of these hidden 
vortices experimentally. 

We use the state shown in Figs. l(a)-l(b) as an exam- 
ple. After a short expansion time, the state begins to look 
very differently. In Fig. 2, the density distribution and 
phase distribution at the expansion time r = 4 is plotted. 
We see clearly in Figs. 2(a)-2(b) that, besides eighteen 
vortices already shown in the in situ density distribution 
[see Fig. 1(a)], a series of new ordinary vortices appear 
along the symmetric axis of two condensates. These new 
visible vortices originate physically from the hidden vor- 
tices. It is not difficult to understand the revelation of 
the hidden vortices during the free expansion. The core 
of a vortex (hidden, ghost, or visible) is also a velocity 
singularity, where the velocity approaches infinity. Be- 
cause the kinetic energy should be finite, during the free 
expansion where the angular momentum is conserved, no 
atoms will be allowed into the core area. As a result, the 
core is stable and will not be destroyed. At the same 
time, as the two BECs begin to overlap, atoms begin to 
move into the central region and fill the space between 
hidden vortex cores, rendering them visible. No all hid- 
den vortices can be revealed during free expansion. As 
shown in Fig. 2(c), there still exist several remnant hid- 
den vortices along the boundary between the two fiows 
[see Fig. 2(c)], which is due to the strong repulsion and 
pushing of the newly formed visible vortices during the 
expansion. Note that in our simulation, we did not find 
the generation of the vortex- antivortex pairs predicted 
in the numerical study of the interference of sliced BECs 
without rotation (ssj . 

Because the ghost vortices lie on the outskirts, the loca- 
tions of their phase defects move outward during the free 
expansion. Consequently, the density in the regime of 
the ghost vortices is always negligible. This is the reason 
that the ghost vortices will not become visible vortices 
during the expansion. 



III. VORTEX FORMATION PROCESS AND 
CRITICAL ROTATING FREQUENCY FOR A 
ROTATING DW BEC 

The vortex formation process with the rotating DW 
potential is drastically different from the one with a sin- 
gle harmonic potential. There is a critical angular fre- 
quency for a rotating single-well potential to create vor- 
tices: when the angular frequency fl is below the crit- 
ical angular frequency, only ghost vortices are formed 
at the outskirts of the BEC cloud and the cloud does 
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FIG. 2: (color online) (a) Density distribution, (b) density 
contour, and (c) phase distribution after the cloud free ex- 
pands for T — 4. Before the free expansion, the system has 
rotated with SI — 0.9 for t = 250. The darker color area or 
contour indicate the lower density or phase. 



not carry any significant angular momentum; when n is 
larger than the critical frequency, visible vortices begin 
to appear along with a jump in the angular momentum 
[24] . For the rotating DW system, the vortex formation 
starts with a pair of hidden vortices. As seen in Figs. 
3(a)-3(e), the hidden vortex pair begin their formation 
at the ends of the potential barrier, then move towards 
the center. This is followed by a sequence of other hid- 
den vortex pairs. Ghost vortices begin to appear only 
after several pairs of hidden vortices are already formed. 
Eventually at the critical rotating frequency ilc — 0.59, a 
pair of ordinary visible vortices are formed [see Fig. 3(f)], 
along with a jump in the system angular momentum [see 
Fig. 3(g)]. 

The dependence of the angular momentum per atom 
Iz on the rotating frequency 51 is shown in Fig. 3(g). 
It is clear from the figure that the angular momentum 1^ 
increases gradually and continuously with D, until a jump 
occurs at fie — 0.59. Along with Figs. 3(a)-3(e), this 
means that the hidden vortices can gradually increase 
the system angular momentum as they move towards the 
center. The ghost vortex has no capacity to carry angular 
momentum. As demonstrated clearly in Figs. 3(d)-3(e), 
even as a pair of ghost vortices move towards the center 
and eventually become a pair of ordinary visible vortices, 
the change in the angular momentum is very sudden as 
witnessed by the jump in Fig. 3(g). 

In virtue of data fit, another interesting feature in Fig. 
3(g) is that, for fl < Uc, Iz increases linearly; in contrast, 
for > flc, Iz grows exponentially. These two different 
regimes marked by well linear growth and perfectly ex- 
ponential growth are likely associated with the fact that 
the hidden vortices only form along the central barrier 
but the ordinary visible vortices can emerge in the whole 
regions of the cloud. 



IV. HIDDEN VORTICES IN A BEC CONFINED 
IN A ROTATING TOROIDAL TRAP 

Hidden vortex can exist in non-DW potentials. To 
illustrate this point, we consider a BEC confined in a 
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FIG. 3: (color online) Phase distributions of ^ at f = 250 
for different rotating frequencies: (a) Q = 0.1; (b) Q. — 0.12; 
(c) n = 0.5; (d) Q. = 0.58; (e) Q. = 0.59, respectively, (f) 
The density distribution \ipf at f = 250 for Q. = 0.59. (g) 
Angular momentum per atom Iz versus Q. The darker color 
area indicates the lower phase or density. 

rotating toroidal trap. The toroidal trap is given by 

Vu{x, y) = + , (3) 

where a characterizes the anisotropy of the 2D central 
barrier in the toroidal trap, a — 1 corresponds to a 
circular toroidal trap , which has recently been studied 
by Aftalion et aL[3J]. We focus on the anisotropic (de- 
formed) toroidal trap (a 7^ 1), where the lack of the 
rotation symmetry excludes the formation possibility of 
multi-quantized vortex (giant vortex) in the center of the 
trap for sufficiently large rotation frequency VL and suffi- 
ciently narrow barrier. The numerical procedure for this 
toroidal trap is identical to the DW potential. In Fig. 
4 we display the steady density distributions (left) 
and the corresponding phase distributions of i/i (right) at 
t = 250 for the anisotropic toroidal trap rotating with 
VI = 0.6 (top) and Vl = 0.9 (bottom), respectively. The 
parameters are Vq = 40, a = 0.8, cr = 1, 7 — 0.03, and 
c = 600. 

At f2 = 0.6, two visible vortices appear in the in situ 
density distribution as shown in Fig. 4(a). Moreover, 
there is an ellipsoid density hole in the trap center which 
looks like a giant vortex. In the phase distribution dis- 
played in Fig. 4(b), we see that besides two phase sin- 
gularities corresponding to the two visible vortices, there 
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FIG. 4: (color online) Density distributions \^'^ (left) and 
phase distributions of ijj (right) at t = 250 for the toroidal trap 
rotating with Q. = 0.6 (top) and = 0.9 (bottom), respec- 
tively. The corresponding parameters are Vb = 40, a = 0.8, 
o = 1, 7 = 0.03, and c = 600. The value of the phase varies 
continuously from to 2-k. The darker color area indicates 
the lower density or phase. 



are other phase defects, which are distributed along the 
long axis of the central barrier and at the outskirts of the 
cloud. The four singly quantized phase defects along the 
long axis of the central barrier show that the ellipsoid 
density hole is not a giant vortex. Our numerical sim- 
ulation further indicates that the four single-quantized 
phase defects carry angular momentum and satisfy the 
Feynman's rule together with the two visible vortices. 
Therefore, they are four singly quantized hidden vortices. 
On the other hand, the phase defects at the outskirts of 
the cloud are ghost vortices because they contribute no 
angular momentum to the system. 

With the increase of rotation frequency, more vortices 
nucleate and a triangular vortex lattice forms eventually 
[see Fig. 4(c)]. At the same time, as shown in Fig. 4(d), 
more hidden vortices also show up in the central barrier 
region, e.g., there are six hidden vortices for VL — 0.9. 
However, these hidden vortices do not form triangular 
lattice, and they are still distributed along the long axis 
of the ellipsoid central barrier. This shows from another 
perspective that the hidden vortex is different from the 
usual visible vortex. The local enlargements of Figs. 4(b) 
and 4(d) are given in Fig. 5, where the circles denote the 
positions of hidden vortices. 
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FIG. 5: (color online) (a) Local enlargement of Fig. 4(b) and 
(b) that of Fig. 4(d), where the circles mark the positions of 
hidden vortices. 



V. DISCUSSION AND SUMMARY 

In summary, we have investigated numerically the for- 
mation of vortices in a rotating double-well potential. 
We found that, other than usual visible vortex and ghost 
vortex, there exists another type of vortex, which we call 
hidden vortex. Unlike the usual visible vortex, these hid- 
den vortices are invisible in the in situ density distribu- 
tion. They differ also from the ghost vortex by being 
able to carry angular momentum. In addition, the core 
size of hidden vortex is not given by the healing length, 
and is strongly influenced by the shape of external po- 
tential. Only after including the hidden vortices can the 
Feynman's rule be satisfied. 

Hidden vortex has appeared in literature in other 
names. For example, the magnetic fluxons in a super- 
conductin g lo ng Josephson junction in a parallel mag- 
netic field [26| , Josephson vortices between two long par- 
allel coupled atomic BECs [l^l , and rotational fluxons of 
BECs in rotating coplanar double-ring traps [28| . The gi- 
ant vortices (or, sometimes, called "phantom vortices") 
in a cylindrical hard-walled bucket or a quadratic plus 
quartic trap [11] or a circular toroidal trap [s^ can also 
be regarded as a form of hidden vortex. However, as we 
have illustrated in the last section with an anisotropic 
toroidal trap, hidden vortex can occur in many settings 
other than previously mentioned structure or potentials. 
Therefore, hidden vortex is a more general notion that 
encompass all the essential features of Josephson vortex 
and giant vortex. At the same time, these names, such 
as Josephson vortex and giant vortex, each coined for a 
special potential, show that it is necessary to distinguish 
hidden vortex from the usual visible vortex and the ghost 
vortex. 
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